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O ■ Abstract 
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We discuss the problem of gauge invariance of the vector meson photoproduc- 
tion at small x within the two-gluon exchange model. It is found that the gauge 



133 invariance is fulfilled if one includes the graphs with higher Fock states in the me- 

son wave function. The obtained results are used to estimate the amplitudes with 
, longitudinal and transverse photon and vector meson polarization. 

(N ' 

^ \ Investigation of vector meson photoproduction at small x is a problem of consider- 

able interest. We are interested in the low- x region where the predominant contribu- 
■ tion is determined by the two-gluon exchange and the vector meson is produced via the 

photon-two-gluon fusion. The factorization of diffractive vector meson production with 
longitudinally polarized photons into the hard part and parton distribution was shown 
in [1]. Thus, such processes, can be an excellent tool to study the generalized parton 
distribution [2]. Moreover, they should give important information on the vector meson 
wave function. The spin-density matrix elements which were studied at DESY (see [3] and 
references therein) should be sensitive to the vector meson wave function. To analyze spin 
effects in the 7* — > V transition, it is necessary to calculate the amplitude with transverse 

. ^ I polarization of a vector meson. For the light meson production, this transition amplitude 

is not well denned because of the present end-point singularities [4]. One of the possible 
ways to regularize these end-point divergences is to include the transverse quark motion, 
as it was done, e.g., in [5, 6, 7, 8]. 

Unfortunately, such higher-twist effects can result in loss of the gauge invariance (GI) 
of the amplitude. In this report, we study the 7* — > V transition amplitude for different 
polarization of photon and vector meson at small x and check the GI of our results. The 
vector meson production can be described in terms of the kinematic variables which are 
follows: 

q 2 = (L - L'f = -Q\ rl = {P- P'f = t,x P = q ' {P q ~/'\ s = (q + Pf, (1) 

where L, L' and P, P' are the initial and final lepton and proton momenta, respectively, 
Q 2 is the photon virtuality, rp is the momentum carried by the two-gluons, xp is part 
of proton momentum carried by the two-gluon system and s is the photon-proton energy 
squared. The vector meson is produced by the photon-two-gluon fusion and the momen- 
tum V = (q + rp) is on the mass shell. The xp variable which is equivalent to skewedness 
( is determined by 

M*+Q 2 + \t\ 

Xp ~ Q ~ — v - . (2) 
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Within the two-gluon exchange model we calculate the L — > L, T — > T and T — > L 
amplitudes which are of importance in analyses of spin density matrix elements. In 
calculations the k- dependent wave function [9] is used 

* v = g[(f+ M v )$y + J- fM - A( f _ M V ){E V ■ K)}(f) v (k, r). (3) 

Here V is a vector meson momentum and My is its mass, E v is a meson polarization vector 
and K is a quark transverse momentum. The first term in (3) represents the standard wave 
function of the vector meson. The leading twist contribution to the longitudinal vector 
meson polarization is determined by the MyJfiy term in (3). The k- dependent terms of 
the wave function are essential for the transverse amplitude of the light mesons. The wave 
function (3) has quite a general form and can reproduce results of most models [6, 7, 8]. 
The other model for the wave function which has a structure similar to (3) was considered 
in [10]. The GI of the vector meson production amplitude was discussed in [5, 11]. It was 
found that the 7* — > V transition amplitude at zero momentum transfer should vanish as 
l\ for l\ — > 0, where l± is the transverse part of the gluon momentum. The importance 
of the higher Fock states of the wave function in GI of the vector meson production 
was shown in [11]. These results were obtained in the two-gluon model exchange in the 
Feynman gauge. 

The leading term of the amplitude of diffract ive vector meson production is mainly 
imaginary. The imaginary part of the amplitude can be written as an integral over z and 
k±. The leading over s term of the 7* — > V amplitude has the form 

Tx v ^ =N dz dkl ,VV \] ±; , ( 4 ) 

J J (k 2 ± + Q 2 ) (k 2 ± + \t\ + Q*) 

where N is normalization, Q 2 = m 2 q + zzQ 2 , z = 1 — z and m q is a quark mass. Generally, 
the numerator of the hard scattering amplitude A\ v ^ can be written as follows: 

Axy^ = A° Xy Xri + A l Xy X ^ l\. (5) 

Only the second term in (5) obeys the GI and appears in (4). The imaginary part of 
the vector meson production amplitude (4) depends on the generalized gluon distribution 
J^{X = C, ..). It can be connected with the unintegrated gluon distribution G through 
the integration over l± 

r l 2 , <k 2 . +Q 2 + \t\ (PL (I 2 ) 

J?(C,f,*i + Q 2 + |*|)= / G(Zi,C,...)- (6) 

nj Jo (11 + A 2 )((Z ± + f ± ) 2 + A 2 ) V± ; V; 

Here rj_ is the transverse part of the rp momentum, A is some effective gluon mass. The 
distribution J-q(x, 0, 5q) is normalized to (xg(x,q 2 )). The l\ factor in the numerator of 
(6) appears from the second GI term of (5). 

Unfortunately, in the model with the higher twist effects like the transverse quark 
motion, the sum of the graphs where gluons are coupled with the quarks in the loop does 
not obey GI. Let us discuss this problem in detail for the L — > L amplitude. The GI term 
of the amplitude has the form 



A 1 - 4— 



Q 2 + k 2 ± (l-Azz) -2m q M v zz] (Q 2 + k 2 ± ) . (7) 
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For the gauge-dependent term (GDT) we have 

A\ L = 2-^= [kl{l - Azz) + m q K - 2M v zz)\ (Q 2 + k 2 ± f . (8) 
v W 

It can be seen that in the nonrelativistic limit z — z — 1/2, m q = My/2 the GDT A° L L is 
equal to zero. For light quarks, when m q = 0, the A° L L term is equal to zero at k\ = 0. 

At the same time, the A° L L term has additional power of (Q 2 + k 2 ^ that compensates 
one propagator in (4) with respect to the GI term (7). As a result, the GDT A° L L of the 
amplitude (4) is similar to the contribution of the higher Fock state. Really, here one 
gluon is coupled directly to the wave function and one quark propagator disappears. Let 
us suppose that we can write the sum of GDT and the contribution of the higher Fock 
state in the form 

Kl ~ < L + B(z, k\) $^(1 + C iL). (9) 

Here by the C l\/Q 2 term in (9) we estimate the higher twist contributions in the qqg term 
of the wave function. Let us suppose that the higher Fock term B(z, k"j_) Q qqg compensates 
the A\ L term in (9). In this case, the contribution proportional to l\ in A LL can be 
estimated as 

Kl~-C 1 ^A\ l . (10) 

One can see that this GDT will be suppressed with respect to the GI contribution (7) as 
a power of Q 2 . Really, 

Al,l rn 2 q + k 2 



«^7^ (11) 



11 Kl q 2 

and we have GI of the model at sufficiently high Q 2 . 

Similar calculations have been done for the amplitude with transversely polarized 
photons and vector mesons. The GI term of this amplitude has the form 

A% T ~ ^-Q 2 \k\i\ + Azz) + My (2M v zz - m q {\ - Azz))] (e[e\). (12) 

My L - 1 

For light meson production the resulting amplitude is proportional to k\. For heavy 
mesons, the term proportional to M v appears too. In the transverse case, the GDT 
which does not vanish as l\ in (5) takes place like for the longitudinal amplitude. If we 
suppose the same compensation of GDT as in (9), we find 

oc zz. {16) 



l\ A l T T 

This means that in the transverse case we do not find a Q 2 suppression of additional 
GDT, but we have only its numerical suppression. Really, it can be seen that the T t ,t 
amplitude has additional divergence like \j{zz) with respect to the T LL amplitude (4). 
In the At,t GDT the additional zz term in the numerator cancels this divergence and 
leads to the numerical suppression of the GDT contribution. 

The GI term of the T — > L transition amplitude is determined by 
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It can be found that in this case we have the numerical suppression of a possible GDT 
contribution like for the TT amplitude (13). 

Thus, we have found that the GDT in the 7* — > V transition amplitudes are suppressed 
and one can use the GI terms (7), (12) and (14) to calculate spin-dependent amplitudes 
of the vector meson production. The average momentum transfer, which is used in (4, 
14) is about < \t\ >~ 0.13GeV 2 [3]. The corresponding amplitudes were calculated for 
the k— dependent wave function (3) with the exponential form of <f>v(z, k\) [12] 



<f>v(z, k 2 ± ) = H exp (- 



2 U 2 
V 



kjb 



zz 



)■ 



(15) 



Here H is a normalization factor. Transverse momentum integration of (15) leads to 
the asymptotic form of a meson distribution amplitude <&y S = 6zz. The model has one 
parameter by which determines the average value of k\ and provides the regularization of 
the integrals in the end- point region. In our calculation, we use the value by ~ 0.65 GeV^ 1 
which leads to a reasonable description of the <jl cross section for p production [12]. Then 
the average < k\ >~ 0.6 GeV 2 . 
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Figure 1. Q 2 dependence of the ratio of he- 
licity amplitudes |Tn I /[Too I [3]. The full curve 
- our calculation; dashed line -results of model 
[6], short dashed and dotted lines- models [7] 
and [8] respectively. 
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Figure 2. Q 2 dependence of the ratio of he- 
licity amplitudes | T01 1 / 1 Too I extracted from 
HI and ZEUS measurements of the spin den- 
sity matrix elements in [3]. Lines are the 
same as in Figure 1. 



The results of calculations for the ratio of helicity amplitudes \T TT \/\T LL \ are compared 
in Fig. 1 with the data extracted in [3] from HI and ZEUS measurements of the spin 
density matrix elements. It can be seen that experimental results are reproduced by the 
model quite well. The model gives a reasonable description of the ratio R = olIot- The 
results of the models [6, 7, 8] are shown in this graph too. We can see that all the models 
describe experimental data satisfactorily. 

The comparison of model results for the \T LT \/\T LL \ with experiment is presented in 
Fig. 2. It can be seen that the T^ T amplitude is more sensitive to the structure of the wave 
function. The best description of the \Tlt\/\Tll\ ratio is found in our model and in the 
models [6, 8]. Note that the experimental errors in the spin-density matrix elements are 
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quite large. This does not allow us to find out which model of the wave function describes 
experiment data adequately. 

In this report, the results of the model for the 7* — > V transition amplitude which 
considers the transverse quark motion have been analyzed. These higher twist effects 
regularize the end-points singularities of the amplitudes but lead in the models to violation 
of GI. Note that a similar problem with GI should take place in the models [5, 6]. It is 
found that the contribution of GDT in the model should be small. This permits us to use 
the model results for the GI terms of the 7* — > V amplitudes for numerical calculations. 
Our results describe experimental data on the ratio of helicity amplitudes quite well. 
Unfortunately, the experimental errors in DESY experiments are large and all known 
models describe the experimental results qualitatively. To obtain more information on 
the form of the vector meson wave function it is important to reduce the experimental 
errors. We hope that the precise analyses of spin density matrix elements can be done in 
the COMPASS experiment at CERN. 

This work was supported in part by the Russian Foundation for Basic Research, Grant 
03-02-16816. 
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